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ANTISYMMETRY IN THE FARADAY ROTATION SKY CAUSED BY A NEARBY MAGNETIZED BUBBLE 

M. Wolleben 1 ' 12 , A. Fletcher 2 , T. L. Landecker 1 , E. Carretti 4 , J. M. Dickey 5 , B. M. Gaensler 6 , M. 
Haverkorn 3,11 , N. McClure-Griffiths 4 , W. Reich 9 , and A. R. Taylor 10 

ABSTRACT 

Rotation measures of pulsars and extragalactic point sources have been known to reveal large- 
scale antisymmetries in the Faraday rotation sky with respect to the Galactic plane and halo that 
have been interpreted as signatures of the mean magnetic field in the Galactic halo. We describe 
Faraday rotation measurements of the diffuse Galactic polarized radio emission over a large region 
in the northern Galactic hemisphere. Through application of Rotation Measure Synthesis we achieve 
sensitive Faraday rotation maps with high angular resolution, capable of revealing fine-scale structures 
of ~ 1° in the Faraday rotation sky. Our analysis suggests that the observed antisymmetry in the 
Faraday rotation sky at b > 0° is dominated by the magnetic field around a local HI bubble at a 
distance of 100 pc, and not by the magnetic field of the Galactic halo. We derive physical properties 
of the magnetic field of this shell, which we find to be 20 - 34 /j,G strong. It is clear that the diffuse 
polarized radio emission contains important information about the local magneto-ionic medium, which 
cannot yet be derived from Faraday rotation measures of extragalactic sources or pulsars alone. 
Subject headings: magnetic fields — polarization — ISM: bubbles — ISM: magnetic fields — Galaxy: 
halo 



1. INTRODUCTION 

Linearly polarized radio emission is Faraday rotated 
as it passes through the interstellar medium (ISM) by an 
angle 

rnearside 

A9 = 4>\ 2 = 0.81 / n c {l)B\\(l)dl A 2 , (1) 

J farside 

where <fr, called the Faraday depth, is the integral along 
the line-of-sight (l.o.s.) ?/[p c ] °f the product of the free 
electron density n e /[cm~ 3 ] and the component of the 
magnetic field along the l.o.s. , Bn /[ /iG] . For wavelengths 
(A/[m]) greater than a few cm, Ad is a few degrees or 
more for typical interstellar values of n c « 0.1 cm" 3 and 
Bh « 1 //G along kpc pathlengths. Thus, in principle, 
the Faraday rotation of polarized emission can be used 
to identify the structure of the magnetic field of the Milky 
Way. 

The Faraday rotation of (mainly unresolved) polarized 
extragalactic sources (EGS) and pulsars can be readily 
measured but provides an irregularly sampled, sparse 
grid of <fi which may include contributions from Fara- 
day rotation within the sources. In the case of EGS the 
Faraday rotation arising in the Galactic magneto-ionic 
medium is integrated along the entire l.o.s. through the 
Milky Way, while using pulsars one can try to compare 
the Faraday rotation from pulsars that are nearby on 
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the sky but lie at different distances to sample the mag- 
netic field in localized regions. In either case, the exis- 
tence of small-scale, high-RM structures in the Galactic 
neighbourhood makes it difficult to distinguish the con- 
tributions of large local features, such as the North- Polar 
Spur (NPS), from the global Milky Way magnetic field 
structure. The interpretation of the data requires care- 
ful statistical analysis and modelling to try to tease apart 
these effects. 

In the case of polarized synchrotron emission originat- 
ing from the diffuse ISM, the Faraday rotation is uni- 
formly sampled at the scale of the telescope beam but 
a major problem is that Faraday rotation can depolarize 
the emission, with a non-linear dependence on frequency, 
so <fi cannot always be reliably identified with the rotation 
measure 

RM = Af?/ A A 2 , (2) 

obtained from observations at different wavelengths (see 
Burn 1966t lSokoloff et al.l[l99l . This can be partly over- 
come using Rotation Measu re Synthesis (RM-synthesis, 
iBrentjens Sz de Bruvnl[2005l ) with observations that con- 
tain many frequency channels across a wide band. 
Choosing a value for cf>, the observed polarization vector 
in each frequency channel is derotated by ^A 2 and the 
rotated vectors are coherently added to obtain an image 
of polarized intensity at that value of Faraday depth. 
The final product of RM-synthesis is a 3-D data cube (of 
polarized emission) with the first and second dimensions 
being the coordinates on the sky, and the third dimen- 
sion being Faraday depth, cf> [rad m -2 ]. Faraday depth 
is usually not related to physical distance. 

We have produced an RM-synthesis cube of a large 
field toward the Galactic centre, using new wide-band 
data obtained with the 26-m Telescope of the Dominion 
Radio Astrophysical Observatory (DRAO). In this paper 
we present evidence that a nearby HI bubble produces a 
large-scale RM pattern on the sky that mimics the ef- 
fect that an overall antisymmetry in the Milky Way's 



2 



M. Wolleben et al. 



Faraday Rotation Map (Dominant-RM) 




80 60 40 20 -20 -40 -60 



GLON <deg> 

Fig. 1. — A Faraday rotation map showing the Faraday depth <j> 
of the dominant emission in each pixel (the position of the peak in 
the RM-synthesis spectrum for each pixel) . The grey scale is chosen 
to make structures with positive <j> more visible. The contour line 
corresponds to <f> = ra d m — 2 . Four yellow lines, fitted by-eye, 
indicate the location of the four filaments that can be identified in 
this map, and which are located in the eastern part of the system 
of HI shells. 

regular magnetic field would have on EGS and pulsar 
Faraday rotation, both in magnitude and sign. Our re- 
sults demonstrate the power of RM-synthesis and high 
resolution polarization surveys in separating global from 
local magnetic field structures in the Milky Way. 

2. OBSERVATIONS 

The polarization data presented in this paper are the 
first results from the Global Magneto-Ionic Medium Sur- 
vey (GMIMS3). The frequency coverage is 1277 to 1740 
MHz with the band split into 2048 frequency channels. 
The angular resolution at the present sky coverage is 
~ 1°. The rms noise at this stage of the survey is 25 
mK in a single channel, and is 1 mK in an RM-sy n thesis 
frame . The reader is referred to IWolleben et al.l (2009, 
2010) for more information on observing and calibration. 

RM-synthesis was performed on data cubes of Stokes U 
and Q. Incomplete frequency coverage causes sidelobes 
in the Rotation Measure Spread Function (RMSF, the 
response function with which the RM-synthesis spectra 
are convolved) which were partially removed by applying 
a CLEAN algorithm to the data cube. The expected 
width of the RMSF is 132 rad m~ 2 (note that structures 
at different Faraday depths can be distinguished with 
higher accuracy). By using a narrower RMSF for the 
cleaning p rocess a final resolution of 60 rad mT 2 was 
achieved 1 !. 

Intrinsic polarization angles cannot be recovered accu- 
rately from our RM-Synthesis data due to the low S/N 
ratio of the emission associated with the bubble. F or a 
detailed discussion see iBrentjens fc de Bruvril (|2005l ). 

2.1. Errors and uncertainties 

Systematic errors in the data are caused by two uncer- 
tainties in the calibration. First, the aperture efficiency 
of the DRAO 26-m Telescope, which is used to calcu- 
late main beam brightness temperatures in U and Q, is 

11 https : //www. astrosci . ca/users/drao/gmims 

12 We carefully checked that all conclusions derived from the 
super-resolved RM cube also apply to the un-CLEANED RM cube. 
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Fig. 2. — RM-synthesis frames showing polarized intensity at 
Faraday depths of 60 rad m -2 (a), rad m — 2 (b), and -50 rad m — 2 
(c). The four yellow lines from Fig. [TJ are repeated in (a). Blue 
and red lines in (b) and (c) indicate polarized filaments that were 
identified in RM-synthesis frames. 
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Fig. 3. — Map of the HI brightness temperature at ui sr = 
15 krns" 1 in a logarithmic grey scale. The position of the polar- 
ized filaments, as well as the position of the Upper Sco subgroup 
(Sco OB2_2) today ("*") and 5 Myr ago ("+") are shown. The 
brightness temperatures of the HI shells are about 3 K stronger 
than the background. 

known only at 1.4 GHz. For the data calibration we as- 
sumed this value to be constant over the band, neglecting 
possible variation of up to a few percent. The effect of 
a 20% gradient of aperture efficiency over the band on 
the final RM-synthesis frames was tested and, while the 
systematic error in the polarized intensities was found to 
be < 10%, the shape and position of polarized structures 
were not affected on a level relevant for this analysis. 
The second uncertainty is the absolute calibration (zero- 
levels) of U and Q, which was derived by subtracting 
ground emission profiles obtained by averaging over a re- 
gion of the sky relatively free of Galactic emission - a 
method shown to work reasonably well for long scans, if 
confusion with Galactic emission can be avoided. 

In order to estimate the reliability of the calibrated 
data we compared polarized intensity (PI ) with an older, 
preci sely calibrated survey at 1.4 GHz (jWolleben et alJ 
2006). We find a good correlation down to PI « 100 mK. 
Since the correlation can be expected to suffer from band- 
width depolarization in the wide-band data, we use this 
value as a conservative estimate of the reliability thresh- 
old in our RM-synthesis data. In the present analysis we 
consider only polarized structures above that threshold. 

2.2. Filaments of strong polarized emission associated 
with HI shells 

Figure [1] shows a Faraday rotation map obtained from 
the RM-synthesis cube. Shown in this map is <j> of the 
dominant emission in each pixel (i.e. the Faraday depth 
of the strongest polarized emission, which is not necessar- 
ily the highest or lowest <f> along that l.o.s.). The global 
properties of the Faraday rotation shown in this map re- 
semble the large-scale pat tern of pulsar and EGS RMs 
seen in the same field fe.g. lHan et aLlll997t rTavlor et al.l 
2009) in both sign and magnitude of the Faraday ro- 
tation. The high sensitivity to local polarized emission 
as well as the unprecedented angular resolution of our 
data, however, allow us to also clearly identify small-scale 
structures. We see filamentary structures of high posi- 
tive (j) in the region 10° < I < 60° and 0° < b < +50°, 
embedded in a large region of positive RMs. 

Figure [2] shows three RM-synthesis maps at = 



+60 rad m~ 2 , rad m -2 , and —50 rad m -2 . The emis- 
sion in these maps was Faraday rotated by the corre- 
sponding 4> somewhere along the l.o.s. Note that these 
are not the only RM-synthesis frames that contain po- 
larized emission from this region: the frames shown in 
Fig. [2] have been selected in order to highlight filamen- 
tary structures at different Faraday depths. However, 
there is polarized emission with Faraday depths up to 
±100 rad m" 2 . 

The filaments seen in Fig. Q] appear again in the 
+60 rad m~ 2 map as emission structures but also in 
the rad m -2 map as depressions: the </> = rad m 
frame shows a "negative imprint" of the filaments and 
the amplitudes of the polarized emission in the two 
frames in this region is similar. The missing emission 
in the rad m -2 frame has simply been shifted into the 
+60 rad m~ 2 frame by the action of a Faraday screen. 

More polarized emission filaments are identified at 

= rad m~ 2 and —50 rad m -2 . The three filaments in 
the —50 rad m~ 2 map (red lines in Fig. are strongest 
in polarized intensity at ~ —20 rad m~ 2 but are most eas- 
ily distinguished from the background in the map shown 
here. 

Figure [3] shows HI at a velocity of lSkms -1 , obtained 
from the LAB-HI survey (jKalberla et al.ll2005h at a res- 
olution of ~ 0.6° for the same field as Fig. [T]and [2] At 
this velocity, shells of HI emission forming a bubble (or 
bubbles) are clearly visible. Projected onto the sky the 
bubble is elongated in Galactic longitude with a major 
axis of about 70° , minor axis of about 40° , and centred at 

1 = 10° and b = +25°. The filaments of polarized emis- 
sion (Fig. Q] and [2|) closely follow the shape of the HI 
shells, with the connection becoming even clearer when 
scrolling through the HI emission in neighbouring veloc- 
ity channels. 

In the following Sections we discuss the physical prop- 
erties of this 2200 deg 2 HI bubble and deduce its mag- 
netic field structure from our RM-synthesis cube. Note 
that the region we are discussing covers 1/20 of the sky. 

3. DISTANCE TO AND ORIGIN OF THE FILAMENTS 

Large, extra-planar, expanding HI shells in the di- 
rection of the Galactic centre are associated with stel- 
lar winds and supernova remnants (SNR) originating 
in the stars o f the Scorpius-Centaurus OB association 
(de Geus 1992). Of the three identified subgroups in Sco- 
Cen, two have been identified with the o rigin of the Lo- 
cal B ubble and the Loop I superbubble (Maiz-ADellaniz 
[200U> . The third sub- group, Upper Scorpius (also know 
as Sco OB2 2), is the furt hest away at a distance of 
145 pc (|de Zeeuw et al.lli~999l ) and its position of I = -8°, 
b = +21° places it in the centre of the western side of 
the HI shells shown in Fig. [3l 5 Myr ago, when this 
subgroup formed, it w as located at I — 15°, b = +21° 
(|Marz-Apellanizl l200lT) : the first OB stellar winds and 
SNR in the Upper Sco subgroup, expected to occur when 
it was 3-5 Myr old, will have occurred as it travelled 
along the major axis of the system of HI shells. Fur- 
thermore, 5 Myr ago the 14-15 Myr old subgroup Upper 
Centaurus Lupus ( UCL) was located at I = 4°, b = +9° 
about 100 pc away (|Maiz-Apellanizl I2001D . The HI fila- 
ments in the north-east of Fig. [3] could plausibly have 
been swept up by SNR and winds from the UCL sub- 
group. 
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If the centre of the HI bubble is at a distance of 145 pc 
then its linear size is about 200 pc x 100 pc. If the cross 
section along the l.o.s. is circular then the nearside of 
the bubble is 95 pc away and the far side 195 pc distant. 
These dimensions are compatible with other, indepen- 
dent, information that constrains the bubble size: since it 
shows no sign of breaking out of the disk its height above 
the mid-plane must be less than the warm HI scaleheight 
of about 300 pc; the depolarizing effect of the 170 pc dis- 
tant HH-region Sh2-27 (centre at I = 6.1°, b = +23.6°, 
diameter of 11°) can be seen in FigJSJb) so most of the 
polarized emission in this region must originate beyond 
170 pc; the wall of the local cavity is about 80 pc distant 
in this direction so the HI shells of this bubble should be 
no nearer than this. 

Mixed emitting and Faraday-rotating regions produce 
polarized emission from a continuous range of Faraday 
depths because emission from the near edge undergoes 
little Faraday rotation whereas emission from the far side 
is rotated by the whole layer. The observed RM-spectra 
for l.o.s. in this region generally have a strong peak (with 
polarized intensities > 100 mK), which we interpret as 
the Faraday rotation of the background polarized emis- 
sion, but many l.o.s. show weak emission from a range 
of Faraday depth, indicating that low level (< 50 mK) 
emission might also be coming from the shell. The shell 
acts as a Faraday screen to the strong background po- 
larized emission (thin in <j)) and also as a much weaker 
mixed emitting and rotating slab. No radio continuum 
emission from this bubble is seen in total intensity at 
408 MHz (lHaslam et al.| [l982l) and 1.4 GHz (|Reichlll982t 
iReich fe Reichlll986l ). but ~ 50 mK emission is too low 
to stand out from the synchrotron background in these 
maps. 

The Faraday depths of the emission rotated by the 
bubble allow us to determine its magnetic field configu- 
ration. Note that with 1° resolution all structures within 
the bubble that are bigger than about 3 pc are resolved. 
The positive (j> in the east (Fig. [2^) changes to negative <j> 
in the west (Fig. [2J:), suggesting that the magnetic field 
is wrapped around the bubble, pointing towards us in 
the east and away from us in the west. The polarized 
emission at 4> = along the centre of the bubble links 
these two regions: = means that either Bn = or 
n c = 0. In this central region the magnetic field pattern 
indica ted by optic al polarization measurements of local 
stars ([Heiles 2000) , as well as radio polarization data at 
1.4 GHz (| Wolleben et al.ll2006h . show that B ± closely fol- 
lows the shape of the HI bubble. This suggests that, in 
this part of the shell, Bn must be small with a B± com- 
ponent of unknown strength. The bottom of the bubble 
lies near the Galactic plane and is not clearly visible in 
any <j) frame. 

Theoretical models of magnetic fields in interstellar 
bubbl es symmetrically inflated by stellar winds and SNR 
(e.g. lFerriere et al.lll99l1 ; lTomisakalll992D predict that the 
ambient B-field is compressed in a shell by the expand- 
ing shock wave. The shape and implied magnetic field 
configuration of this HI bubble suggest that it has ex- 
panded asymmetrically. In particular, for Bn to have 
constant magnitude but opposite signs at the two ends 
of the bubble indicates that its expansion is constrained 
to only one direction along the l.o.s. The dense wall of the 
Local Bubble lies in the right direction and at the right 
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Fig. 4. — The sketch depicts the position and geometry of the 
HI shell and its magnetic field relative to the Local Bubble. Sco 
OB2_2 is drawn at its current position. GC denotes Galactic cen- 
tre. 

distance to constrain the expansion of shells powered by 
the Upper Scorpius sub-group (5-6 Myr old) towards us. 
It seems likely that magnetic pressure in the wall of the 
Local Bubble has caused an assymmetric expansion. Fig- 
ure 4 shows a sketch of the deduced bubble and magnetic 
field configuration. 

4. DISCUSSION 

We have searched the co mbined 

VTSS/SHASSA/WHAM survey dFinkbeinerl I2003T) 
for signs of the HI shells and bubble in Ha emission 
but could not find any. Assuming no optical extinction, 
this suggests that the high Faraday rotation from the 
shells is caused by a local enhancement of magnetic 
field and not electron density. The emission measure 
towards the shells is around 2.4 cm -6 pc. Assuming a 
scale height for the War m Ionized Me dium of 0.9 kpc 
(iSavaee fc Wakker I [20091 ) to 1.8 kpc (jGaensler et al. I 
[2001 we have an average electron density along b = 20° 
of 0.03-0.02 cm -3 . The shells have an approximate 
thickness of 5° or Sr ~ 13 pc, giving a maximum path 
length through a shell of L &s 2\ / 2rSr ks 100 pc. At the 
ends of the major axis of the bubble, where we expect 
the magnetic field to be mainly directed along the l.o.s., 
the observed Faraday depth of \<j>\ « 50-60 radm -2 cor- 
responds to a magnetic field strength of Bn « 20-34 /iG. 
It is noteworthy that the Faraday rotation pattern 
has been fractured into confined filaments; naively one 
would expect a Faraday rotation distribution which 
changes smoothly from positive on one side to negative 
on the other. 

Such strong magnetic fields are difficult to reconcile 
with a simple sweeping up of a typical ambient 2-3 fiG 
field by a single strong adiabatic shock, such as a SNR 
in the Sedov- Taylor phase, which has a compression ra- 
tio of only four. Two possibilities to explain the strong 
magnetic fields are: (i) the SNR and/or stellar winds 
that inflated the bubble have entered the radiative stage 
of their evolution (plausible given the 5-6 Myr age of 
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the sub-group), allowing stronger compression to oc- 
cur; (ii) the Upper Sco sub-group OB stars had al- 
ready entered the dense wall of the local bubble, with 
its already compress ed ma gnetic field of about 8-9 /zG 
(jAndersson fc Potted 1200ft iMao et~aD I2010f >. when the 
new bubble expanded. The latter possibility is compati- 
ble with the elongated shape (Section 3). 

The main astrophysical consequence of this nearby 
bubble and its associated Faraday rotation is its effect 
on the interpretation of EGS and pulsar Faraday rota- 
tion in terms of the global magnetic field structure of the 
Milky Way. 

Maps of pulsar RMs (e.g. lHan et al.llT997f ) show a clear 
antisymmetry in the sign of RM in the northern Galactic 
hemisphere looking toward the Galactic centre, with pos- 
itive RMs dominating in the region 60° > I > 0°, b > 0° 
and negative RMs at 0° > I > -60°, b > The recent 
Farad ay rotation map of NVSS sources by iTavlor et al.l 
( 2009), which provides the most densely sampled grid of 
RMs to date, shows a similar antisymmetric pattern of 
RMs. The RMs of pulsars and EGS in this region have 
magnitudes up to 100radm~ 2 , similar to the maximum 
Faraday depths found in our RM-synthesis cube. 

It is tempting to interpret such large-scale an- 
tisymmetries in Faraday rotation as a signature 
of the regular magnet ic field o f the Mil ky Way 
(lAndreasyan fc Makarovj 11988k lHan et alj 119971 : 
ISun et al.l I2008T ). especially as galactic dynamo theory 
makes predictions about the ex pected mean field sym- 
metries in the disk and halo (e.g. lBeck et al. Ill99ft) . But 
several authors have pointed out that the antisymmetry 
may be due to radio L oop I and its bright, polarized 
segme nt, the NPS. iSimard-Normandinfc Kronberd 
(1980) find that the effect of the NPS on the RMs in 
the northern Ga l actic hemisphere is clearly visible and 
llnoue fc Tabaral (jT981) point out the importance of 
avoiding the effect of irregularities in RMs due to local 
objects. 

In this paper we have identified a nearby bubble and 
set of shells that covers a large area of the sky with the 
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same Faraday rotation asymmetry as seen in the EGS 
and pulsar RMs. The discovery of small-scale structures 
in the Faraday rotation sky which are associated with 
nearby HI, have the same sign, similar magnitude, and 
are embedded in the large-scale antisymmetry at b > 0° 
leads us to conclude that the observed antisymmetry has 
largely a local origin. Note that the structures discussed 
here are not associated with the NPS as they do not re- 
semble the shape of the spur: thus the NPS can only add 
to the region of the northern sky where local structures 
dominate the Faraday rotation pattern. To minimize the 
effect of local RM structures on the picture of the mean 
Galactic magnetic field we suggest that point source RMs 
towards regions that show diffuse polarized emission at 
high Faraday depth should not be used. 

Clearly the diffuse polarized radio emission contains 
important information about the local magneto-ionic 
medium, with better angular resolution than currently 
available from EGS source measurements. This must be 
taken into account when building a complete picture of 
the Galactic magnetic field. At present RM-synthesis 
data exist only for the northern sky. Observations for 
a southern sister-survey (STAPS, principal investigator: 
M. Haverkorn) were just completed and RM-Synth esis 
should soon become available too. Ide Geusl (|1992h (in 
their Fig. lb) show more HI filaments associated with 
Sco OB2_2 in the southern sky. Our prediction is that 
there will be Faraday rotation structures also associated 
with these filaments. 
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